Opaque or transparent? A link between neutrino optical depths 
and the characteristic duration of short Gamma Ray Bursts 
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ABSTRACT 

Cosmological gamma ray bursts (GRBs) are thought to occur from violent hypercritical ac- 
cretion onto stellar mass black holes, either following core collapse in massive stars or compact 
binary mergers. This dichotomy may be reflected in the two classes of bursts having different 
durations. Dynamical calculations of the evolution of these systems are essential if one is to 
establish characteristic, relevant timescales. We show here for the first time the result of dynam- 
ical simulations, lasting approximately one second, of post-merger accretion disks around black 
holes, using a realistic equation of state and considering neutrino emission processes. We find 
that the inclusion of neutrino optical depth effects produces important qualitative temporal and 
spatial transitions in the evolution and structure of the disk, which may directly reflect upon the 
duration and variability of short GRBs. 
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The coalescence of two compact objects due to 
the emission of gravitational radiation is a pos- 
sible origin of cosmological gamma ray bursts of 
the short variety (see e.g., Meszaros 2002, for a re- 
view), lasting a few tenths of a second (Lattimer 
& Schramm 1976; Paczyhski 1986; Eichler et al. 
1989; Narayan, Paczyhski & Piran 1992). Pre- 
vious multidimensional studies have shown that 
the outcome of such a merger is in all likelihood 
a dense torus surrounding a supramassive neu- 
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tron star, likely to collapse to a black hole on a 
short timescale, or a black hole, if one was already 
present in the system (e.g Rasio & Shapiro 1994; 
Kluzniak & Lee 1998; Ruffert & Janka 1999; Ross- 
wog, Ramirez-Ruiz & Davies 2003). This class of 
systems has been studied observationally since the 
1970s, when the Hulse-Taylor system was discov- 
ered (Hulse & Taylor 1975), and the recent detec- 
tion of PSR J0737-3039 (Burgay et al. 2003), the 
tightest binary yet in this class, has only height- 
ened interest in these events. The final instants 
before the actual merger will produce a powerful 
burst of gravitational waves, and the merger wave- 
form itself is expected to reveal details about the 
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equation of state of matter at high densities (see 
e.g., Thorne 1995). 

The merger process is intrinsically a three- 
dimensional event, with no particular symmetry 
which can be exploited to reduce the complexity 
of the problem, and must be addressed numeri- 
cally. No such calculation is currently evolved for 
more than a few tens of milliseconds, because of 
computational limitations. Information about the 
state of the system after the violent, dynamical 
merger is essentially over can be obtained, but the 
final state of the system, on a timescale compa- 
rable to the duration of a short GRB, cannot be 
addressed. 

The steady state structure and composition of 
post-merger accretion disks has been the focus of 
several recent studies (Popham, Woosley & Fryer 
1999; Narayan, Piran & Kumar 2001; Kohri & Mi- 
neshige 2002; Di Mattco, Pcrna & Narayan 2002; 
Pruet, Woosley & Hoffman 2003; Beloborodov 

2003) , which have included an increasing amount 
of microphysical detail. However, the dynamics of 
the problem have never been followed for an ex- 
tended amount of time, although progress is being 
made in this direction (Setiawan, Ruffert & Janka 

2004) . 

In a previous paper (Lee & Ramirez-Ruiz 
2002), we studied the dynamical evolution of the 
accretion structures formed after a black hole- 
neutron star merger under a set of simplifying as- 
sumptions. Initial conditions were taken directly 
from the late stages of the coalescence, modeled 
in three dimensions using a simple equation of 
state (Lee 2001). The disk was then mapped onto 
two dimensions assuming azimuthal symmetry, 
and its evolution was followed for 0.2 seconds, 
comparable to the duration of a short GRB. The 
equation of state was simple (ideal gas with adi- 
abatic index 7 = 4/3), and we assumed that all 
the energy dissipated by viscosity (modeled with 
the full expression of the stress tensor, using an 
a-prescription for the magnitude of the viscosity) 
was radiated away. 

In this Letter we present for the first time 
dynamical calculations of post-merger accretion 
disks in two dimensions, (r, z), that: (1) use a re- 
alistic and self-consistent equation of state for the 
fluid in the disk; (2) consider the cooling of the 
disk through neutrino emission, using appropriate 
rates; (3) account for the effects of neutrino opaci- 



ties and the corresponding optical depths through 
a simplified treatment and (4) last for a minimum 
of 0.4 seconds, comparable to the duration of a 
short GRB. 

2. Physics 

After the initial merger phase of the binary is 
completed, a few tenths of a solar mass can be 
left in orbit around the central black hole (which 
has a mass Mbh ~ 3-5 M Q ), in a small disk a few 
hundred kilometers across. The densities are high, 
with 10 9 < p(g cm -3 ) < 10 12 , and the correspond- 
ing temperatures are 10 10 < T(K) < 10 11 . Under 
these conditions, nuclei are entirely photodisinte- 
grated, and the fluid consists of free electrons, pro- 
tons, nuclcons, plus a small number of positrons. 
The nucleons make up a non-relativistic ideal gas, 
while electrons are degenerate and highly rela- 
tivistic. Photons arc completely trapped due to 
the extremely high optical depth, and are hence 
advected with the flow. The composition of the 
fluid is determined self-consistently assuming: (1) 
charge neutrality and (2) complete /3-equilibrium 
(Kohri & Mineshige 2002; Lee, Ramirez-Ruiz & 
Page 2004). This allows us to determine the elec- 
tron fraction Y e and the electron degeneracy pa- 
rameter (through the chemical potential) locally, 
and to follow it throughout the dynamical evolu- 
tion. This is important since the composition di- 
rectly affects the neutrino emission rates (see be- 
low). The equation of state is then given by 

3 m u VMe/ 

where K is the constant corresponding to a fully 
degenerate, relativistic electron gas, /i e = 1/Y e 
is the number of baryons per electron, and P v is 
the pressure due to the presence of neutrinos in 
the optically thick regions of the disk (see below). 
Strictly speaking, the electrons arc not fully de- 
generate, since rj e /kT is in the range 2-4 (rj e is 
the electron chemical potential). However, this 
represents a small error in the equation of state, 
since gas pressure dominates the overall balance, 
contributing about 90% of the total pressure, fol- 
lowed by electron degeneracy pressure at around 
5%. 

Neutrinos are emitted abundantly in the disk, 
due to the high temperatures and densities. We 
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take into account all of the following processes: (1) 
nucleon-nucleon bremmsstrahlung; (2) pair an- 
nihilation; (3) electron and positron captures onto 
free nucleons and (4) plasmon decay. In practice, 
the global cooling rate is completely dominated 
by captures, and is hence well approximated 
by (Kohri & Mineshige 2002): 



Qcap 



1.1 x 10 3 



kT 



T{ x erg s" 1 cm" 3 (2) 



in the regime where electrons are degenerate. This 
expression has been used directly to compute the 
cooling within the disk. The error in doing so, as 
opposed to a full integration over the phase space, 
is limited to about a factor of 2 (in the inner re- 
gions, where the cooling is most intense), as com- 
pared with rates extrapolated from the work of 
Langanke & Martfnez-Pinedo (2001). 

The emitted neutrinos arc not entirely free to 
leave the system, since scattering off free nucle- 
ons is important. The optical depth for neutrinos 
can be estimated as H/l v , where H is a typical 
disk scale height, and l v is the mean free path 
between scatterings. The cross section depends 
on the mean neutrino energy as a n cx E% (Tubbs 
& Schramm 1975), and E v « 43(r e/ 9i 2 ) 1/3 MeV 
(P12 = p/10 12 g cm -3 ), using the fact that neu- 
trinos are being produced by reactions involving 
degenerate electrons (see e.g., Shapiro & Teukol- 
sky 1983). This amounts to a surface of last 
scattering for neutrinos, or "neutrino-surface" at 
p ~ 10 11 g cm~ 3 . Our simple treatment of this 
fact consists of suppressing the local cooling rate 
by a factor exp(— t„), and adding a pressure term 
P v oc aT 4 [l — cxp(— t,,)] to the equation of state. 
This simple alteration is crucial in determining the 
energy output of the disk. We note that no modifi- 
cation in the composition was computed owing to 
this fact. In reality, the opaqueness of the material 
will lead to an increase in the electron fraction Y e , 
by up to a factor of 2 (Beloborodov 2003). The 
optical depth will thus rise by a factor 2 2 / 3 ~ 1.6, 
making the optically thick region in the disk larger 
and the overall luminosity slightly lower. Any en- 
ergy release in neutrinos will thus be spread over 
a longer period, making our estimates concerning 
time scales fall on the conservative side (see below, 
Section 3). 

The dynamical evolution is followed for at least 
0.4 seconds in the central Newtonian potential 



produced by the black hole, $ = —GMbu/R, us- 
ing a two dimensional smooth particle hydrody- 
namics (SPH) code in azimuthal symmetry (Mon- 
aghan 1992). Accretion is modeled by placing an 
absorbing boundary at the Schwarszchild radius 
ra = 2GMbh/c 2 , and the mass of the black hole is 
updated continuously. There is no external agent 
feeding the disk with matter, and no boundary 
conditions need to be specified (other than for ac- 
cretion). The equations of motion and the energy 
equation contain all the terms for a physical vis- 
cosity derived from the stress tensor (Flebbe et al. 
1994), and we use an a-prescription for the coef- 
ficient of viscosity, with 10~ 3 < a < 10 _1 . The 
dissipated energy is injected back into the disk as 
thermal energy, and may leave the system depend- 
ing on the local cooling rate (see above). 

3. Temporal and spatial transitions 

If the disk is assumed to be optically thin every- 
where, and the effects discussed above are ignored, 
the neutrino luminosity is a smooth, monotoni- 
cally decreasing function of time, similar to what 
we obtained in previous work (we have confirmed 
this in trial runs with the present input physics). 
Accounting for a variable optical depth leads to 
behavior that is markedly different, as can be seen 
in Figure 1, where L v (t) is shown for three values 
of the viscosity parameter a. For a high viscos- 
ity, a = 0.1, the accretion timescale is so short 
(40 ms), that most of the material is accreted 
onto the black hole and has no time to radiate 
away its internal energy reservoir of « 10 52 erg. 
The disk becomes entirely optically thin at around 
30 ms, and this explains the rcbrightening that is 
observed. The luminosity subsequently exhibits a 
power law decay, as in the trial runs with no opti- 
cal depth effects, with L v oc i~ 9 . 

For lower viscosities (a = 10~ 2 , 10~ 3 ), the be- 
havior is qualitatively different at early times, and 
is a consequence of the longer accretion timescales: 
0.4 s and 1 s respectively (we note here that all 
our estimates of t acc are based on the actual ac- 
cretion rate, i.e. t acc *=s M<j/Mbh)- Quite simply 
the material in the disk is allowed to remain in 
the vicinity of the black hole without being ac- 
creted, at high temperatures and densities (since 
the transport of angular momentum is so ineffi- 
cient), and radiate away essentially all of its in- 
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ternal energy through the emission of neutrinos. 
Once this occurs, on a cooling timescale given by 
tcooi ~ Ei n t/L u , a break occurs (at t ~ 80, 100 ms 
for a = 10~ 2 , 10 -3 respectively) and the luminos- 
ity drops again as t~ 9 . At this stage the en- 
ergy output comes directly from viscous dissipa- 
tion within the disk. The break time is determined 
by how much energy the disk holds initially, E- lnt , 
and by the rate at which it is lost, L v {t), which 
is itself fixed by the densities and temperatures 
found within it. The limitation on the cooling 
rate imposed by high optical depths is essential, 
and allows the energy loss to be spread over an 
extended period of time. For the three calcula- 
tions shown here, the total energy release in neu- 
trinos, up to 0.4 s, is {E u /10 52 crg)= 1.2, 1.05, 0.8, 
for a = 10 _1 , 10~ 2 , 10~ 3 , respectively. Assuming 
that this energy is converted to e ± with a one 
per cent efficiency in the regions along the rota- 
tion axis through vv annihilation, approximately 
10 50 erg could be released and made available for 
the production of a relativistic fireball. A detailed 
calculation of this is left as future work. 

The other important mechanism by which a 
GRB might be powered is MHD energy extraction, 
either through the Blandford-Znajek mechanism 
(Blandford & Znajek 1977), or through a magnet- 
ically dominated jet. Since our code does not ex- 
plicitly account for the presence of magnetic fields, 
we make a simple, rough estimate of the magnetic 
field strength by assuming that its energy density, 
B 2 /8ir, is in equiparition with the internal energy 
density, pc 2 (this is clearly an upper limit) . Within 
the disk, the field will be most intense where the 
density is highest. Figure 2 shows the maximum 
density in the disk as a function of time, for the 
same calculations as shown in Figure 1. The in- 
ferred magnetic field is B ~ 10 15 - 10 16 G. It 
is again evident that the high-viscosity case is 
markedly different from the other two, showing a 
significant prompt decrease in /9 max . The transi- 
tion to complete transparency to neutrinos occurs 
at around 40 ms. The balance between cooling, 
which produces vertical compression (and hence 
a rise in central density), and accretion, which 
slowly reduces the mass of the disk, produces a 
break in the curves for a — 0.01,0.001 at t « 
0.05, 0.8 s respectively. The total estimated energy 
release through the Blandford-Znajck mechanism 
is (£ B z/10 51 erg)= .34,1.3,6.5, for a = 10" 1 (up 
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Fig. 1 . — Neutrino luminosity as a function of time 
for three runs using identical initial conditions and 
different values of the viscosity parameter: a = 
10" 1 , a = 10~ 2 and a = 10~ 3 . 
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Fig. 2. — Maximum density as a function of time 
for the same runs as shown in Figure 1 . The inter- 
val marked AT S corresponds to the sound crossing 
time across the optically thick region of the disk. 
The subsequent rise is due to the vertical compres- 
sion of the disk as it cools. Accretion eventually 
drains the disk of matter, producing the character- 
istic breaks at t = 0.05,0.8 s for a = 10~ 2 , 10~ 3 . 



4 



to 0.4 s), a = 10~ 2 (up to 0.4 s), a = 10~ 3 (up to 
1.2 s), respectively. 

In addition to the temporal variations in total 
luminosity, there is a spatial transition within the 
disk, when the optical depth is around unity. For 
the calculations with a < 10~ 2 , there is always 
(up to the end of the run) an opaque region in 
the center of the disk. The density p and entropy 
per baryon s/k along the equatorial plane (z = 0), 
as well as density contours in the r — z plane, are 
plotted in Figure 3, 0.2 s after the beginning of the 
calculation with a = 10~ 2 . We also show for com- 
parison the equatorial density profile of the disk 
in a calculation that did not take into account the 
effects of optical depths (dashed line in Figure 3). 
There is a clear break at r « 9 x 10 6 cm, where 
t v = 1. As the optical depth increases and the 
cooling is suppressed, the pressure rises, inhibit- 
ing the thinning of the disk and the increase in 
density in the inner regions. This break is essen- 
tially a stationary feature, due to the relatively 
long accretion time scales, as long as the optically 
thick region is present. We believe this transi- 
tion is significant, as geometrically thick accretion 
flows around rotating black holes, with H ~ R, are 
believed to be more favorable for the production of 
energetic outflows than thin flows, where H <C R 
(Livio, Ogilvie & Pringle 1999; Meier 2001). 

4. Discussion 

The evolution of accretion disks resulting from 
dynamical three dimensional binary coalescence 
calculations, where a neutron star is tidally dis- 
rupted before being swallowed by its black hole 
companion, is studied numerically. Angular mo- 
mentum transport and the associated energy dis- 
sipation are modeled using an a prescription. By 
assuming azimuthal symmetry we are able to fol- 
low the time dependence of the disk structure for 
a fraction of a second, a time comparable with the 
duration of short GRBs. During this time, the 
rate of mass supply to the central black hole is 
of the order of a fraction of a solar mass per sec- 
ond; i.e. much greater than the Eddington rate. 
Although the gas photon opacities are large, the 
disk becomes sufficiently dense and hot to cool via 
neutrino emission. There is in principle no diffi- 
culty in dissipating the disk internal energy, but 
the problem is in allowing these neutrinos to es- 
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Fig. 3. — Density profile and entropy per baryon 
along the equatorial plane, z — 0, 0.2 s after the 
start of the calculation with a = 10 -2 (solid line). 
The dashed line indicates the density profile for 
the corresponding trial run that did not include 
optical depth effects. The vertical mark on each 
curve indicates the position where r v = 1. In- 
set Panel: Logarithmic density contours (equally 
spaced every dex). The neutrino-surface corre- 
sponds roughly to logp = 10 11 g cm~ 3 . Under 
these circumstances, the gas cannot efficiently cool 
and mantains a pressure scale height comparable 
with its radius. This, however, cannot happen too 
close to the symmetry axis where the centrifugal 
force cannot oppose the vertical pull of gravity and 
a funnel is formed. 

cape from the Thomson thick inflowing gas. 

At sufficiently low accretion rates, a < 0.01, we 
find that the energy released by viscous dissipation 
is almost completely radiated away on a timcscale 
given by t coo i ~ E- m t/L v ^0.1 s. In contrast, for 
a higher mass supply, a > 0.1, energy advection 
remains important until the entire disk becomes 
optically thin, at t[r v = 1] ~ 30 ms. The restric- 
tion on the cooling rate imposed by high optical 
depths is key because it allows the energy loss to 
be spread over an extended period of time (i.e. 
t coo \ or t[r v = 1]) during which the neutrino lumi- 
nosity stays roughly constant. In principle, neutri- 
nos could give rise to a relativistic pair dominated 
wind if they converted into e^- pairs in a region 
of low baryon density. This gives a characteris- 
tic timescale for energy extraction, and may be 
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essential for determining the duration of neutrino- 
driven, short GRBs. 

An alternative way to tap the torus energy is 
via magnetic fields threading the disk: the energy 
liberated by accretion is converted efficiently into 
magnetic form and emitted as a magnetically dom- 
inated outflow. The launching of a jet probably 
requires the existence of a poloidal magnetic field 
of magnitude B p over a scale of radius R, where 
Bp is smaller than the field strength £?disk asso- 
ciated with the dynamo-driven magnetic disk vis- 
cosity (Livio, Ogilvie & Pringle 1999). This re- 
sults directly from estimating the accretion and 
jet luminosities as L acc - M acc t^, and L jet ~ 
BpR 2 v<f, respectively, which gives (Bp/B^isk) 2 ~ 
(Lj ct /L acc )(H/R). Assuming the scaling B p ~ 
(H/R)B disk , derived by Tout & Pringle (1996), 
one obtains Lj et /L a cc ~ H/R. In order to achieve 
Lj ct ~ iacc, Livio, Pringle & King (2003) ar- 
gue that a number R/H of neighboring annuli in 
the disk need to provide a locally generated net 
poloidal field in the same direction. They notice 
that since the local dynamos vary on timescales of 
idyn, the timescale for establishing a change in the 
magnetic field in the disk should be of the order 
£jct ~ t dra 2 R / H . In our case R/H ~ 3, which gives 
tjct ~ 10 — 100 ms. This corresponds to the typ- 
ical variability timescale of short duration GRBs 
(Ramirez-Ruiz & Fenimore 2000; Nakar & Piran 

2002) . 

Two further physical effects, noted previously 
in a more general context, are relevant for the 
present study. First, the mechanism by which a 
jet is launched is not scale-free (Pringle 1993), 
and thus the poloidal field must be established in 
the inner regions of the disk. Second, in these 
same regions, the inflow speed caused by angu- 
lar momentum loss to the outflow is increased, 
up to a factor R/H over the value expected from 
magnetic viscosity alone (Livio, Pringle & King 

2003) . Thus the poloidal flux generated by the 
dynamo can be effectively trapped. This proba- 
bly implies that once the dynamo process gener- 
ates a global poloidal field, it is able to maintain 
it during the lifetime of the thick portion of the 
disk, t[r v — 1]. For the typical parameters stud- 
ied here, this timescale could be of the order of 
seconds when the rate of mass supply is low (i.e. 
a ~ 10~ 3 ). There is thus no problem in principle 
in accounting for sporadic large-amplitude vari- 



ability on timescales as short as one millisecond, 
even in the most long-lived short GRBs. 
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